Kondo Temperature for the Two-Channel Kondo Models of Tunneling Centers 
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The possibility for a two-channel Kondo {2CK) non Fermi liquid state to appear in a metal as 
a result of the interaction between electrons and movable structural defects is revisited. As usual, 
the defect is modeled by a heavy particle moving in an almost symmetric double-well potential 
(DWP). Taking into account only the two lowest states in DWP is known to lead to a Kondo- 
like Hamiltonian with rather low Kondo temperature, Tk. We prove that, in contrast to previous 
believes, the contribution of higher excited states in DWP does not enhance Tk- On the contrary, 
Tk is reduced by three orders of magnitude as compared with the two-level model: the prefactor in 
Tk is determined by the spacing between the second and the third levels in DWP rather than by 
the electron Fermi energy. Moreover, Tk, turns out to be parametrically smaller than the splitting 
between the two lowest levels. Therefore, there is no microscopic model of movable defects which 
may justify non-Fermi liquid 2CK phenomenology. 



It is well known that two-level systems (TLS) deter- 
mine the low energy phenomena in a glassy matter. The 
most popular realization of the TLS is a movable atom 
tunneling between two minima of the two-well poten- 
tial created by other atoms The low-temperature 
behavior of glasses was found to be consistent with the 
assumption of homogeneous distribution of both energy 
difference and spatial distance a between the minima. 
In metallic glasses TLS interact with itinerant electrons. 
Usually this interaction in metallic glasses is assumed to 
be weak and to manifest itself only in a finite relaxation 
rate of the TLS, see Rcf. [01 for a review. 



It was proposed long ago that a TLS interacting 

with itinerant electrons behaves like a localized spin in 
the Kondo model. Indeed, in the limit kpa <ti 1, where 
kp is the Fermi wavelength, only the electrons with two 
spherical harmonics, namely, I = and I — l,m — 
interact with TLS. (Here and below the axis of the mo- 
mentum quantization is the easy axis of TLS, x) 

Let us introduce a pseudospin of a symmetric TLS: 
S — —1/2 corresponds to the ground state (even wave 
fimction), whereas S = 1/2 labels the excited state 
with the odd wave function. One can map the elec- 
trons with the two relevant spherical harmonics on the 
one-dimensional Fermi gas of particles that are charac- 
terized by a pseudo-spin with components cr = -|-, — as 
^1=0 = ^'/=i,m=o = ^+7 while the real electron spin 
index s =t, i is replaced with the channel index /i — 0, 1. 
Furthermore, provided that the Fermi energy, ep, ex- 
ceeds all of the relevant energies, one can linearize the 
electron dispersion law near the Fermi level, e{p) ~ i'fIpL 
where vp is the Fermi velocity. The divergences caused 
by the linearized spectrum should, thus, be cut off by 
the bandwidth D ep. The resulting Hamiltonian of 
the system can be expressed as 
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Here the Pauli matrices, t^^, , act in the space of the elec- 
tron pseudospin, and S is the operator of the TLS pseu- 
dospin, [S\S^] = ie'^''Sk. The first term in Eq. (0) de- 
scribes kinetic energies of ID electrons. The second term 
characterizes the TLS level splitting: and repre- 
sent correspondingly the tunneling and the initial TLS 
asymmetry. The third term in Eq. (^) describes TLS- 
electron interaction. The Hamiltonian Eq. is nothing 
but the two-channel Kondo Hamiltonian M, where the 
level splitting plays the role of the Zeeman splitting of 
states of the usual Kondo impurity. For a comprehensive 
review of implication of this model to magnetic ions and 
tunneling centers in metals see Ref. 

The two-channel Kondo effect is known to manifest it- 
self through a non-Fermi liquid behavior of the specific 
heat, magnetization and electronic correlation functions. 
Such a behavior takes place when both the temperature, 
T, and the level splitting. 
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do not exceed the Kondo temperature, Tk- 
shown M3M that in the limit <C ^ 1 
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The non-Fermi liquid behavior of the TLS at the two- 
channel Kondo fixed point was used in Ref. ||l^] to in- 
terpret the zero bias anomaly in characteristics of point 
contacts |jl^ and more recently |jl^ to explain the tem- 
perature behavior of the dephasing rate observed in 
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Refs. [HQ. In general, TLS were assumed to play im- 
portant role in crystalline metals as well as in metallic 
glasses. 

Assumptions of Ref. [|l2| were criticized on the grounds 
that the disorder induced splitting estimates were too 
high for the Kondo-like behavior to develop Al- 
though several questions raised in Refs. [ p5[ remain unan- 
swered, we put this issue aside. Instead, we concentrate 
on a different objection - smallness of Tk in Eq. (||) for 
reasonable values of parameters. Here we prove that 
Tk A2 for any set of the microscopic parameters, 
which allows the Kondo like description, Eq. (|^). There- 
fore, the two channel Kondo fixed point by no means can 
be reached with the lowering T and thus is irrelevant for 
the description of the TLS in metals. 

To understand why the resulting Kondo temperature 
is so small, let us first discuss the physical meaning of 
the bare coupling constants Vj and estimate them. 

The coupling constant Vx in Eq. determines the 
renormalization of the TLS asymmetry by the electrons 
- a tilting of the double-well potential by the dipole mo- 
ment of the electron density. Assuming a contact interac- 
tion characterized by a dimensionless coupling constant 
A < 1 we can estimate Vx at given Fermi wave number, 
kp, and the size of TLS, a, as (see, e.g., (2]J^) Vx — Xkpa. 

As to Vz, it characterizes transition between the two 
states of the TLS assisted by an electron transition. 
Incoming electron renormalizes the barrier's height, V ^ 
and consequently the tunneling amplitude. However, 
the tunneling event still has to occur. Therefore, [^|J|] 
Vz — Xk^pO? ex\i{—ri), where the tunneling exponent 77 is 
determined by V and the atomic mass M . Since M is 
large, 77 3> 1 even for relatively low barriers. As a result, 
the coupling constant Vz is always much smaller than '■ 

Vz — Vx{kFa)exp{—7]) <t:. Vx , 77 ~ h~^aVSMV (4) 

This is why Vz was usually neglected in previous treat- 
ments of TLS. For "typical" values of the parameters 

D^5eV, -0.2, Vz/vx-lQ-\ (5) 

the "conventional" estimate of the Kondo temperature is 
~ 10~2 - 10-3 K 0]. This low value of Tk makes 
it hard to believe that the Kondo fixed point is relevant 
for the discussion of existing experiments. 

In attempt to resolve the problem of small Tk authors 
of Ref. Q went beyond the two level approximation and 
considered virtual tunneling through the third level of the 
"TLS" . This contribution to Vz apparently does not con- 
tain the tunneling exponent. According to Ref. this 
fact dramatically increases Tk comparing to the "con- 
ventional" estimate. In our opinion the statement about 
the large increase in the Kondo temperature is incorrect. 
Below we discuss the problem in detail. 

Tk can be extracted from second order perturbation 
theory in the interaction of the tunneling particle with 



the electrons. We calculate this correction within the 
one-dimensional (ID) model suggested in Ref. 

Consider a heavy particle in a symmetric ID double- 
well potential V{x). Let the energies of the two lowest 
eigenstates, i?i,2, be indistinguishable: E2 — Ei 0. 
One can express the matrix element of the contact inter- 
action of this particle with electrons through the coupling 
constant, A, Fermi velocity, vp, the heavy-particle eigen- 
functions, 4>i{x), and the electron wave functions, ipa-{x), 
with a given isospin, cr = +, — : 

U:^""' = 27tvf\I-/^' , III' = J dx(j>kM.i^p ■ (6) 

There are two second order corrections to the scatter- 
ing amplitude <7i,i — s- (72, j, which correspond to pro- 
cesses with different intermediate states: 

(i) In the intermediate state an electron has an isospin 
a and the particle occupies a state k with the energy 
Ek = Ei+ek (see Fig. §1). 




FIG. 1. Second order diagrams for the scattering ampli- 
tude, a) electron process, b) hole process. 

(ii) The transition i — > fc of the heavy particle produces 
an electron-hole pair - an electron in the final state (T2 and 
a hole in the state a. Afterwards the hole annihilates the 
electron in the initial state ci and the heavy particle 
changes its state from k to the final state j (Fig. ^3) 

Combining the contributions of these two processes 
and taking into account the occupation numbers of elec- 
tron states at a given temperature T we present the sec- 
ond order correction to the matrix element as: 

^^g -°'^'-SL (-e:-0[l + exp(-,/Tr 

Here I^' is determined by Eq. (||), ^ denotes the energy of 
the intermediate electron or hole state counted from the 
Fermi level, and Ek = Ek — Ei is the energy of the fc-th 
state in the double- well potential counted from its ground 
state El — E2- We have already discussed that the do- 
main of the integration over ^ should be |^| < D ~ £p. 
As to the singularity at ^ = — et, the integral in Eq. (|^) 
should be understood as the principal value. The minus 
sign in front of the second (hole) term in the numerator 
is due to the anti-commutation of fermionic operators. 

Consider now matrix elements 11^2'^'^ = —1/21'^^ that 
describe transitions of the heavy particle between its two 
lowest states. Using Eqs. (|^) one can show that such a 
transition should be accompanied by the change of the 
electron isospin: due to the parities of the wave functions 
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= U^2^ = 0. To evaluate 11^2^ we sum over a and 
integrate over ^ in Eqs. (^. The result can be written as 

n 2 
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where e^^ = e^^ = 0, e^^ = — e^^ = 1, and 



(9) 



We start with the contribution of the first two levels 
(k = 1,2) to the sum over k in Eq. (||). Using Eqs. (|^), 
and (^) one obtains, cf. with Ref. 



U2 = In (^^^ 1 dx [<Pl - 0?] _ ^2 ] 



(10) 



The electron wave functions are standing waves, 

ip-{x) + i'ip+{x) = \/2exp(ifcFa:). Since kpa <C 1 (other- 
wise this ID approach is not applicable) 

^-{x) V2[l- ikFx)^/2] , ip+{x) V2kFX . (11) 

We introduce the wave functions, (/'i(r), localized in the 
left (right) well, \/20i,r — 0i ± 02, recall that the func- 
tions 0i are normalized, and rewrite Eq. (nOf) as 



U2 = 81n^^^ j dx {kpxf (l)i{x)(j)r{x) . (12) 

The wave functions (j^iix) and (j)r{x) are localized in the 
different wells, their overlap being exponentially small. 
Accordingly, U2 is exponentially small as well, and one 
arrives at the "conventional" estimate for Tk- 

It turns out that taking into account the higher excited 
states in the double- well potential, i.e., terms with k > 2 
in the sum Eq. (^), can only reduce the estimation of Tr-. 

Indeed, it follows from the definition of 
Eq. (||), and completeness of the set of the functions, 
{^feji {Y.k't^k{x)(j)k{y) = 5(x - y)) that the sums, 
Sfc ^ik'^" Ijk'^^ ■^ith any set cr„ and consequently J2k 4j' 
are symmetric with respect to the permutation of i and 
j, i-e-, Efe 4 = Sfc 4^- As a result 
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The sum rule, Eq. (|13|), together with Eq. ([g) imply that 
the second order correction Eq. (||) to the matrix element 
Ui2^ vanishes at high temperatures, when — T. More 
precisely, the usual Kondo logarithmic temperature de- 
pendence of Ui2^ persists only for T < = E3 — Ei. 
Therefore, in the expression for Tk, Eq. (|3|), the band- 
width D ^ ep ~5eV should be substituted by £3 ~3meV: 



Tk = £3 {v^v,)^^^ ( ^ 



(14) 



It means that the Kondo temperature is about three or- 
ders of magnitude less than the "conventional" estimate! 

To interpret this result note that when it tunnels, the 
heavy particle is under the barrier for a time ^ h/e^. At 
energies bigger than £3 wc thus deal with a continuously 
moving particle rather than a pseudospin. 

However, any truncation of the sum, Eq. (^), results in 
a strong overestimation of the Kondo temperature. In- 
deed, the truncated sum Wn<oo, Eq- (|l^), is neither zero, 
nor exponentially small. Thus u„ = w„ ln(£'/r) 3> U2 at 
least when £„ <^ T. It is the substitution of Mqo by M3, 
proposed in Ref. Q, that dramatically enhanced Tk- 

The sum rule Eq. ( p^ ) hints that although contribution 
of each excited state k to u„ is quite large in absolute 
value (if k is not too big), these contributions have dif- 
ferent signs and cancel each other up to an exponentially 
small quantity Uqo when all of them are included. 

To demonstrate that this is the case we repeated nu- 
merical calculations of Ref. using the same model 
potential. Fig. ||, but took into account all of the excited 
states {k} rather than only A; = 3. 

Following Ref. |^] we chose the barrier height to be V = 
9.86?i^/2M&2, where M is the particle mass and b is the 
well width. We computed the eigenfunctions (jnix) and 
used Eqs. ( ^),( ^ and (^) to evaluate u„ = Un(kpb)^'^ . 
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FIG. 2. Symmetric double well potential with the well 
width 6, the barrier has a hight V and a width 2ab. 



Fig. ^ shows n-dependence of the ratios j/„ = Unjui = 
Un/u2 (for £f = 10^£3, T = 0.00204£3, and the relative 
width of the barrier a equal to 2.5). One can see that 
2/3 3> 1, i.e., U3 ^ U2. As we expected, absolute values of 
2/4, ?/5, 2/6, are also large, but the signs alternate. In agree- 
ment with our analytical conclusions further increase of n 
gradually reduces |2/n|, and 2/„ — )■ 2/00 ~ 1 when n ^ 00. 

On the insets of Fig. ^ we present a-dependencies of 
2/„ to make it evident that although m„>2 is not ex- 
ponentially small as U2 is, it regains this smallness as 
n — !■ 00. Indeed, 2/34 is almost a constant in the interval 
2.5< a <3.0 whereas 2/5 increases with a by factor ~5 in 
the same interval, and ln(2/5) is a linear function of a. 
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FIG. 3. j/n = u„/u2 as a function of the level number n and 
of the relative width of the barrier, a (insets). Upper inset: 
1/34(0). The variation in the interval 2.5< a <3.0 is less than 
10% . Lower inset: linear dependence of ln(y5) on a. Note 
thaty5(3)/y5(2.5)^5. 

We also computed the temperature dependence of 
'U30 ~ Moo for several Fermi energies. This dependence 
is presented on Fig. ^ in scmilogarithmic scale. All four 
curves coincide , i.e., iioo does not depend on sp- More- 
over, as it was expected, the logarithmic dependence, 
Moo oc ln(r) + const, persists only as long as T < £3 
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FIG. 4. The dependence of the quantity u-m = U3o{kFb) ^ 
on ln(r/e3) for different ef shown in the legend. 

One can see that the numerical simulations unambigu- 
ously support the analytical conclusions. Taking into 
account the excited levels does not remove exponential 
smallness of the second order correction to the scatter- 
ing amplitude. A similar problem with the similar solu- 
tion - contribution of the continuous spectrum above the 
barrier to the a-decay rate - is described in the book . 

Returning to the Kondo temperature, Eq. (|lj), we find 
that Tk < 10"^ for the "typical parameters" Eq. (||) and 
an optimistic estimate £3 « 50K. Therefore the Kondo 
model based on movable structural defects is hardly able 



to explain the experiments, |11 l^JT^ . 

Moreover, the estimate (^4|) excludes the very possibil- 
ity of the development of the strong coupling two-channel 
Kondo regime at arbitrary low temperatures: it implies 
that the splitting of the two lowest levels of a TLS , A, 
Eq. (^), always exceeds Tk- Indeed, A > Aj ~ £36^'' ^ 
i;2£3 ~ A(/cFa)^e~''£3. Using the fact that the model is 
applicable only in the limit ^ <C 1 we obtain 
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The same conclusion can be reached for any double-well 
potential model. Therefore a movable defect weakly cou- 
pled with electrons is unable to demonstrate the two chan- 
nel Kondo non-Fermi liquid behavior. Jl8[ . 
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slightly changes index in Eq. ([l5[): 7^7 = l/{Avx) — 
1/2 — Vx- One can see that 7 « 7 ^ 1, since Vx <^ 1. 
Therefore the inequality ( |l5| ) remains valid even after 
taking into account the renormalization of Az. 
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